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Abstract 

Dynamic spin susceptibility is calculated for the t — J model in the para- 
magnetic phase by applying the memory function method in terms of the 
Hubbard operators. A self-consistent system of equations for the memory 
function is obtained within the mode coupling approximation. Both itinerant 
hole excitations and localized spin fluctuations give contributions to the mem- 
ory function. Spin dynamics have a diffusive character in the hydrodynamic 
limit; spin-wave-like excitations are regained in the high-frequency region. 

I. INTRODUCTION 

While the reference antiferromagnetic (AFM) insulator compounds of high-T^ supercon- 
ductors (HTSC's) are well understood in terms of 2-dimensional (2D) isotropic Heisenberg 
model the nature of anomalous spin-dynamics in the doped samples still requires proper 
understanding . 

One of the simplest models invoked to describe the HTSC's is the t — J model, which 
contains the essential physics of Cu02 planes in the superconducting cuprates. The t — J 
(or its extension t — t' — J) model is the low energy effective model obtained from the 
Hubbard model by projecting out doubly occupied sites. As a result the t — J model is 
formulated in terms of the so called Hubbard operators (HO's), which are neither Fermi nor 
Bose operators. This particularity makes it difficult to treat the t — J model within the 
conventional field-theory methods. 
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The various approaches, e.g., slave-boson or slave- fermion methods and dia- 
grammatic technique for HO's []^, have been used to study the spin dynamics within the t—J 
model. In the slave-field approaches the local constraint is usually replaced by the global 
one restricting the validity of the approach. Whereas in any approximation formulated in 
terms of HO's the constraint of no doubly occupancy can be rigorously preserved. 

Recently in Ref. the diagrammatic technique for HO's has been used to calculate 
the spin susceptibility within the Larkin equation |^ for the t — J model. However, in 
the particular case J = (corresponding to t/ — cxo limit of the Hubbard model) the 
contribution from the irreducible part in the denominator of Larkin equation (see Eq. (9) of 
Ref. P]) vanishes, which indicates that the so called kinematical interaction is not properly 
taken into account. 

In the present paper we study the dynamic spin susceptibility for the t—J model by 
applying memory function technique |jl6[ in term of HO's which has been applied recently 



for calculation of the optical conductivity for this model |T0|. We show that there exist 
two different in nature contributions to the memory function. The first one is due to the 
kinematical interaction and comes from the particle-hole excitations in the itinerant hole 
subsystem. While the second one comes from the localized spin fluctuations due to the 
Heisenberg interaction. The existence of these two contributions explicitly shows that there 
is a competition between itinerant and localized magnetism as it has been pointed out in 
Refs. and observed experimentally (see, e. g. Ref. pj]). 

It is found that the low energy {uj 0) spin dynamics has a diffusive character. While 
in the high energy limit (cu — oo ) the spin-wave-like excitations are regained. The mean- 
filed-like (MFL) expression obtained earlier by Kondo and Yamaji's (KY's) theory 
is recovered in this limit. 

The paper is organized as follows. In the next section we formulate the t — J model 
in terms of HO's. In Sec. 3 the general formalism of the memory function approach is 
presented and within the mode coupling approximation the memory function is calculated. 
The self-consistency of the presented approach is discussed in Sec. 4. Sec. 5 summarizes 



our main results. In the Appendix the expression for static spin susceptibihty is derived. 



II. THE MODEL 

The Hamiltonian of the t — J model reads 

h=h,+hj=y: u.xfxf + -hi^^^, - \^^^^\■ (1) 

The first term in the right hand-side of Eq. (|l]) describes an electron hopping between the 
nearest (tjj = t) and next-nearest itij = t') neighbor lattice sites. The second term describes 
the exchange interaction of localized spins Sj between the nearest neighbor sites {Jij = J). 
The HO's in Eq. (|l|) are defined as 

Xf = \^,a){^,P\ (2) 

for three possible states at the lattice site i 

\i,a) = \i,0), \i,a), (3) 

for an empty site and for a singly occupied site by electron with spin cr/2 (a = ±1). In the 
t — J model only singly occupied sites are retained and the completeness relation for the 
HO's reads as 

^r + E^r = i- (4) 

The spin and density operators in Eq. (0) are expressed by HO's as 

S^ = Xf, St = W^X1% n, = Y.Xr- (5) 
where a = —a. The HO's obey the following multiplication rules 

Xfxf = 5,,5^,Xf (6) 

and commutation relations 
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^ = 6,, {6f,,Xf ± Ss^xf) . (7) 

In Eq.(0) the upper sign stands for the case when both HO's are Fermi-hke ones (as, e. g., 
X^'^). The spin and density operators (j^) are Bose-hke and for them the lower sign in Eq.(^ 
should be taken. The HO's are neither Fermi nor Bose operators, the are projected operators. 
These unconventional commutation relations (|^ makes impossible to treat the model within 
the conventional diagrammatic technique. To use the latter one needs to introduce slave 
particles with the constraint of no doubly occupancy. While in the treatment within the 
HO's the constraint is automatically fulfilled. Therefore we treat the problem in terms of 
HO's and use the memory function formalism to determine the dynamic spin susceptibility. 



III. MEMORY FUNCTION 



A. General Formalism 



The dynamic spin susceptibility is defined as 

where Rjj = Rj — Rj and {{A\B))^ denotes the Fourier transformed two-time retarded 
commutator Green function (GF) |T^jr5[ 



oo 

iujt 



{{A\B))^ = -z dte^^\[A{t),B]), (9) 







where Imu > 0, A(t) = exp{iHt)Aexp{—iHt), and (AB) denotes the equilibrium statistical 
average. 

In the paramagnetic state with zero sublattice magnetization an average of the commuta- 
tor {[S^, Sj']) = 25ij{Sf) equals zero. Since just this quantity enters as an initial condition 
{t = 0) in the equation of motion for the GF (H), it is more convenient to construct the 
self-consistent equation for the Kubo-Mori relaxation function (see, e.g., Ref. ||16|| ) 



{{A\B))^ = -i dte'^\A{t),B), (10) 

^ 



where {A(t),B) is the Kubo-Mori scalar product defined as 

{A{t),B) = ['^ d\{A{t-iX)B), (11) 

where P = 1/T (/?, = fc^ = 1). 

The GFs (|) and (p!oD are coupled by the equation 

u{{A\B))^ = {{A\B))^-{{A\B))^=o. (12) 
There are also following useful relations which can be obtained from the definitions (|9|)- (pAD 

iitA\B))^ = iiA\ - tB))^ = {{A\B)U (13) 

izA,B) = {A,-zB) = {[A,B]), (14) 
{A,B) = -{{A\B))^=o, (15) 

where iA = idA/dt = [A,H]. 

By using the above formulas for the dynamic spin susceptibility we obtain 

Xq(w) = Xq-W'^'g(^), (16) 

where Xq = Xq(0) is the static spin susceptibility and $q(co') = ((•S'q ))a;. To calculate the 
spin-spin relaxation function $q(c(;) it is convenient to employ the memory function approach 
of Mori (see, e.g. Ref. [|16|)- We define the memory function M(q, cj) by the equation 

To calculate the memory function we use the equation of motion for the relaxation 
function ([T0| ) 



u;{{S;\S-))^ = X^ + {{^S^\S-)U (18) 
and similarly ((^•S'q ))^ obey the following equation of motion 
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u;{{tS+\S-))^ = {zS^, S-) + {{zS^\ - zS-))^. (19) 

From Eq.(|14D we have {iS^,S~) = {[S^,S~]) = 2/ y/N {3^)6^0 which is zero in the param- 
agnetic phase, that results 

u;%{u) = xq + ^((^^+1 - tS^)h. (20) 
By introducing the zero order GF $q(c<j) = Xq,/^ we rewrite Eq. (^) as follows 

= + (21) 

where we have introduced the scattering matrix 

T^{u) = -{{zS^\-zS^)l-. (22) 



By comparing (|22D to the definition of the memory function ( p!7|) we get the following relation 



between the scattering matrix and the memory function 

T,M = %^ + %^<MW. (23) 

Aq Aq 

A formal solution of the Eq. ( ^3]) by iteration shows that the quantity M(q, u;)/xq(u;) is 
just the irreducible part of the scattering matrix which has no parts connected by the single 
zero order GF ^^{u): 

M{q,u;) = xlT-^iu;) = ((z5+| - ^5"))^ (24) 
Finally, the dynamic spin susceptibility in terms of the memory function can be written 



as 



Xq = -Xq ^77 vf— • (25) 



B. Mode-coupling approximation 



First we express the memory function in terms of the irreducible current-current time- 
dependent correlation function by using the spectral representation for the GF 
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M(q,.) = ge-^^.4£..' ^,j':;; iye-'-V!it)\J,r, (26) 

where the current operator in the site representation is defined as Jj = iSj'. Current operator 
can be written as a sum of two terms 

Jj = Jj + J' = i^t^ Ht] + Hj] (27) 

In Eq.(|27|) the first term comes from the so called kinematical interaction which is due 
to the unconventional commutation relation for the HQ's operators (|^). This term is pro- 
portional to the hopping integral and reads as 

Jj = -Yl ^jrni.Xf^X'^ - X+^X""). (28) 

m 

The second term in Eq.(^) comes from the exchange interaction between localized spins 
and has the form 

= 2 ^ Jjm{SjS^ — S^Sj'). (29) 

m 

To calculate the irreducible time-dependent correlation function in the right hand -side 
of Eq. (p6D we employ the mode-coupling approximation in terms of an independent propa- 
gation of the dressed particle-hole and spin fiuctuations (see, e. g., Gotze et al., [l^). This 
scheme is essentially equivalent to the self-consistent Born approximation in which the ver- 
tex corrections are neglected. The proposed scheme is defined by the following decoupling 
of the time-dependent correlation functions: 

(^^°(t)^r(t)^r^r) ^ (^„;°w^r)(^r(i)^r), (3o) 

{Sat)S-{t)S^S^) ^ {Sat)S^){S-it)S^). (31) 

The cross-correlations like {Ji\{J'^)j') are ignored within the proposed approximation and 
they will be omitted. 

By using the above defined decoupling scheme (|5y) and ( PT| ) and the spectral represen- 
tation for the GF, the memory function (|26|) can be written as 
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M(q, u) = Mt(q, to) + Mj{q, u), (32) 
where Mt{c{,uj) is the contribution from the itinerant hole subsystem and reads as 

Mtici, = ^ E 4c // du^du' [ni., - - ""ti^l'f^",'/^ ^ (33) 

— oo 

where n{uj) = (e^'^ + tkq = tk - tk-q with tk = ^(^Tq + i^'7q), ^ = 4, 7q = 

1/2 [cos g^: + cosgj^] and 7q = cosg^.' cosg^ for 2D square lattice (the lattice constant is taken 
to be unity) and A]^{uj) = —l/7clm.{{X'^X^))^ is a hole particle spectral function which is 
spin independent in the paramagnetic phase. 

The second contribution Mj(q, cj) in Eq.(^2|) comes from the localized spin subsystem 
and is given by 

OO 

TT'^A' ^ ^JJ UJ [UJ — UJ + IT]) 

where N{uj) = {e^'^ — 1)^^ and J^q = Jk — Jl^-^, and Jq = zJ7q. In obtaining (|3^ relation 
(('S'q|'S'q))aj = l/2((5'q whlch Is Valid in the rotationally invariant system has been 

used. 

The real, ReM(q, ti;), and imaginary, ImM(q, c<j), parts of the memory function are odd 
and even functions of u , respectively, and they are coupled by the dispersion relation 

ReM(q,o;) = ird^- ^"^^^^'"^^ (35) 

Therefore, only imaginary part of the memory function should be evaluated. 

C. Asymptotic behavior of Xq('^) 

Now we examine asymptotic behavior of the dynamic spin susceptibility. First we con- 
sider the hydrodynamic limit q ^ and uj ^ 0. In this limit, ReM(q, cj) being an odd 
function of uj vanishes while ImM(q, u;) remains finite. By using Eqs.([33|) and (0) we can 
express it as ImM(q, uo) ~ — Dq^ with D = Df + Dj where 



8 



= T7Xl(^^ktk)^ lim lim / duJi\n'{uJi)\A]^{uJi)A]^^ci{uJi - uj), 

—oo 

oo 

2 r 
= — - V(qVk^k)^ lim lim / duJi\N'{uJi)\lraXk{uJi)lraXk-ci{uJi - u), (36) 

VriV , cj-^Oq-»0 J 



where q = q/g, Vk = dt^/dk, n'{uj) = dn{uj)/duj, and N'{lj) = dN{uj)/duj. Finally in the 
hydrodynamic limit the dynamic spin-susceptibility can be expressed in the usual form (see 
Ref. PI) as 

Xq(^) = Xo ^ (37) 

where D = D/xo is the spin diffusion coefficient and xo is the static uniform susceptibility. 

Unlike to the hydrodynamic limit, in the high energy limit, — »• oo the dominant 
contribution to the memory functions comes from the real part: M(q, uj) ~ m^juj where 
is the first non- vanishing moment in 1/a; expansion of the memory function defined as 



oo 

-i I rfu;ImM(q,cu) = ([iS+,S+]) (38) 

— oo 

Thus in the high energy limit dynamic susceptibility takes the form 



Xq(t^) = -Xq 2 2 (^^) 



where = mq/xq. Let us note, that Eq. ( pUD with expressions of Xq and mq derived in 
the Appendix [see Eqs. ( [A 161 ) and ( |A17| )] reproduces the result for the spin susceptibility 
obtained for the t — J model in Ref. by KY's theory |]T3| which is essentially self- 
consistent MFL approximation. 



IV. SELF-CONSISTENCY OF THE PROBLEM 

The equations (p5D,(p3D, and (Q) are the self-consistent integral equations for dynamic 
spin susceptibility which is obtained by using only mode coupling approximation. These 
equations should be solved numerically by iteration procedure. The static spin susceptibility 
at each iteration step should be calculated by Eq. (|A4|) with M(q, cj) and Xq(^) from 
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the preceding iteration. However some ansatz for Xq(^) ^is the starting point of iteration 
procedure should be defined. Moreover we need to know the hole spectral function entering 
into Eq. (R3[) for the memory function. 



According to well known results for hole spectral function obtained within t— J model |T8 
(a;) can be modeled as 

A^iu) = Z^6{uj + /i - £k) + ^r(c^), (40) 

where is the quasiparticle weight for the excitations with the dispersion e-^ in a narrow 
band of the order J. The second part A'^^'^{uj) is due to the diffusive motion of holes in a 
broad band with bandwidth 2W (of order 8t for 2D square lattice). We model it as 

Ar = N,^cO{W-\uj + fi\), (41) 

where Nine is the density of state for the incoherent continuum and it is coupled to Zk by 
the sum rule 

I E d^M^) = iX^' + Xn = 1 - ^. (42) 

By using Eq.(^) the contribution from the hole coherent motion to the imaginary part of 
the memory function can be expressed as 

ImMr^(q, 00) = t^,ZkZk-,^^fc^^^^^^^5(ek - .k-, - u;) (43) 

k 

To evaluate the second term in the memory function (|3^ as the starting point for Xq('^) 
one can use the MFL expression (|39|), with mq and cjq defined by ( |A8|) and (|A16|) , respec- 
tively. As a result we obtain 

ImMj(q,^) = |^^Bkq{Pkq(^) + Pkq(-^)} (44) 

k 

where 

B., = (45) 

t^k^^k-q 
10 



is an effective vertex function and 

[ UJ UJ J 

(46) 

Eqs.(^3D and Eq.(^^ can be considered as the first iteration for the memory function. 
The hole parameters entering into expressions ( |A8| ) and ( |A16| ) for mq and tUq can be calcu- 
lated from the hole spectral function (|40|) . Whereas xi and X2 defined by Eqs. ( |A10D and 
( |A18|) should be evaluated self-consistently from the dynamic spin susceptibility (|39D . By 
using the above expressions (^31) and (^) for ImM(q, uj) and the dispersion relation (|35|) the 
dynamic spin susceptibility within the first iteration can be calculated from Eq. (|25|) . 



Using the obtained results one can evaluate the spin fluctuation part of the memory (|34D 
and the static spin susceptibility ([A4|) -(|A^) for the next iteration procedure. The iteration 
procedure should be continued until the convergency will be reached. 

V. SUMMARY 



To summarize, based on the t — J model and memory function approach we have derived 
a general representation for dynamic spin-susceptibility (^) in terms of the memory function 
(p6|). Our approach is formulated in terms of HO's and therefore the constraint of no doubly 
occupancy is rigorously preserved. The memory function is calculated by using the equation 



of motion method for two-time retarded GF's |Tj] within the mode coupling approximation 



( p^ - |3^) . The two contributions to the memory function is obtained. The first one (pSf ) 
comes from the itinerant hole subsystem and is due to the kinematical interaction. The 
second one (|3^) comes from the localized interacting spin subsystem. In the limit of small 
concentration of doped holes the latter one gives the main contribution which describes spin 
dynamics characteristic for the Heisenberg model. Whereas in the opposite limit of large 
hole concentration particle-hole excitations characteristic to the itinerant magnetism give 
the main contribution to spin dynamics. We have shown that in the paramagnetic phase 
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there are two regimes in the spin dynamics. In the hydrodynamic hmit (g — > 0, — > 0) 
the spin susceptibihty (^) describes diffusion spin dynamics with the diffusion coefficient 
(p6|), which has essentially two contributions. While in the high-frequency limit {uj oo) 
spin- wave- like excitations described by Eq. (^) are observed. Their dispersion, Eq. ( [A16D , 
obtained in the mode coupling like approximation for the equal time correlation function 
( |A14|) recovers the earlier MFL result [T^ obtained within the KY's [|T3| theory. 



To compare our results with that obtained by diagrammatic methods we would like to 
point out that our approach, based on the general representation for the spin susceptibility 
(p5|), is equivalent to summation of infinite series of diagrams generated by the memory 
function (pB)). The latter one, being calculated in the mode coupling approximation ,Eqs. 
(0) and ([31|), can be schematically represented by two loop-diagrams: the first one of order 
due to the particle-hole loop and the second one of order due to the spin fluctuation loop. 
In Ref. all the contributions in the denominator of the Larkin equation are proportional to 
J and therefore disappears in the limit J = {U —>■ oo). While in our approach contribution 
due to the first loop (for which the kinematical interaction is responsible) remains. Whereas, 
in Ref. the spin fluctuation contribution given by our second loop is neglected while several 
other diagrams beyond our simple one-loop diagram due to the kinematical interaction are 
taken into account. 

At present time it is difficult to justify any of discussed scheme, including our mode 
coupling approximation. To check the validity of our approximation one has to solve nu- 
merically self-consistent equations and compare the obtained results with the experimental 
data. This will be done in a forthcoming publication. 
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APPENDIX: 



In this Appendix we evaluate the static spin susceptibility Xq following Tserkovnikov [jlQ . 



For this purpose, by using Eqs.([l7D, (p!9|), and ( pOD we can rewrite the memory function or 
the irreducible part of the current-current correlation function in the following way proposed 



by Tserkovnikov |T5[: 



{{^S^\ - = i{^S^\ - ^S-))^ - i{tS^\S-)U{S^\S-)Vi{S^\ - iS-)\. (Al) 

Likewise we define the irreducible part of the force-force correlation function as 

and consequently the irreducible part of equal time force-force correlation function can be 
written as 

{^'S^,^'S-)- = {fS^,z'S-) - {^'S^, S-){S^, S-y^S^z'S-). (A3) 
By using Eq. ( |A3| ) and identities (p!3|)-(|r5|) the static spin susceptibility can be expressed 



as 



, - <l'^-^.l>' (A4) 

The first term in the denominator can be written as the third moment of the dynamic spin 
susceptibility 

1 °° 

i^S^,t'S^) = {[^'Sl -tS-]) = - I dco co'lmx^ico) (A5) 

— oo 

and the latter one is equal to the second non- vanishing moment of the memory function 

oo 

{i^S+i^S-y''^ = -- I c/cuu;2lmM(q,u;). (A6) 



-oo 
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The expression ( [A4| ) is an exact representation for the static spin susceptibihty. 

However, to derive an approximate expression for Xq we start from the following identity 

{[^S^,Sq]) = {^'S^,S^). (A7) 

We evaluate the left hand-side of Eq. ( |A7| ) by using the commutation relation for HO's, 
that results 



{[^S^\S^]) = 4zJ(l - 7q){^^i + ^n[\^ - xi}, (A8) 



with the following notations 



q q 



and 



To calculate the correlation function in the right hand-side of Eq. (|A7|) we employ the 
decoupling scheme which is essentially equivalent to the mode coupling approximation but 
for the equal time correlation function. Due to the unconventional commutation relations 
for HO's it is more convenient to use the site representation 



StS^)=Y.e-'''''-i^'S^,Sn, (All) 



where second derivative of 5*^^ reads 



with 



TT _ -r+O -t-0+ I Y+OfyOOj_ y \vO- 

ntjn — [y^j + )y^n ! 

p OZ OZ 0+ OZ OZ 0+ 

ijn 'Jj '-'n '-'i '-'n'-'j 1 
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In obtaining ( |A12| ) we have neglected terms proportional to tJ since they give no contribution 
within the adopted approximation. 

In the sum (|A12| ) only two site indices can be equal. We extract those terms and by 
using the multiplication rules (|^) replace the product of two HO's with the same site indices 
by one operator. On rearranging, in the sum there are no products of operators having the 
same site indices. Therefore in all products operators can be interchanged. (Of course in 
the case of two Fermi operators one has to change the sign of the product). Further, we 
substitute the properly rearranged right hand-site of Eq. (|A12| ) into ( |A11| ) and make the 
following decoupling 

(xr5+,5r)^(xr)(5;,5r) {i^j) 

{SfS-S^, Sf) ^ Sf) + {S^Sr){St, Sf) {z ^ j ^ n) (A14) 

In the above defined decoupling scheme the operators on the same lattice site is never 
decoupled. Therefore within the adopted approximation the local correlations are retained. 

In the momentum space the above defined decoupling scheme results in the following 
expression: 

{est,, S-) ^ S^) = ulx^, (A15) 

where 

ul = U^z'\xi\{l - 7q)[H- A + CAq + 7q], (A16) 
with the following notations: 

A = ^{x2 + ^\Xi\+a-v}, C=^[^, (A17) 

X2 = ^j:<iS^S^)^ (A18) 

q 

^ (t'y tt' 

« = ^]yj^E7q^q, «'= ^^E(7q)'^q, ^=^^E7q7>q, (Al9) 
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and 

Tl 

n^ = {l---n^), n={Xr). (A20) 

Therefore, from Eq. ( [A7| ) by using and ( |A15| ) for the static spin susceptibihty we 
obtain the following representation 



nrtr. 



X, = ^^^^^^ = 3- (A21) 



q q 



Essentially, the equation for static susceptibility ( |A21| ) with expressions for mq ( [A8|) and 



cUq ( |A16| ) coincides with that one obtained in Ref. [T^ and can be evaluated self-consistently 
from the one-particle GF. 
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